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Abstract. Recently, fractional Chern insulators (FCls), also called fractional 
quantum anomalous Hall (FQAH) states, have been theoretically established in lattice 
systems with topological flat bands. These systems exhibit similar fractionalization 
phenomena as the conventional fractional quantum Hall (FQH) systems. Using the 
mapping relationship between the FQH states and the FCl/FQAH states, we construct 
the many-body wave functions of the fermionic FCl/FQAH states on the disk geometry 
with the aid of the generalized Pauli principle (GPP) and the Jack polynomials. 
Compared with the ground state by exact diagonalization method, the wave-function 
overlap is higher than 0.97 even when the Hilbert space dimension is as large as 3 x 10®. 
We also use the GPP and the Jack polynomials to construct edge excitations for the 
ferminoic FCl/FQAH states. The quasi-degeneracy sequences of fermionic FCl/FQAH 
systems reproduce the prediction of the chiral Luttinger liquid theory, complementing 
the exact diagonalization results with larger lattice sizes and more particles. 
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1. Introduction 

The integer quantum Hall (IQH) effect [I] and fractional quantum Hall (FQH) effect |2] 
have broken the conventional framework of condensed matter theories, and have dehned 
new topological states of matter. In 1988, Haldane [3] proposed a prototype lattice 
version of the IQHE without an external magnetic field, and thus without Landau levels 
(LLs), which has two non-trivial topological bands labeled by Chern numbers |1]. The 
topological states in Haldane-model-like systems are also called as quantum anomalous 
Hall (QAH) states or Chern insulators (CIs). Based on the understanding of the 
CI/QAH states, and considering the strongly correlated fermionic or bosonic systems, 
possible FQAH states or fractional Chern insulators (FCIs) can be further conceptually 
defined O |6l [TJ El El [101 HB 113]. In the absence of an external magnetic field, the 
topological flat bands (TFBs) [HI El [7] in CI/QAH systems play the role of LLs, and 
can also host the intriguing fractionalization phenomena. Since the TFBs quench the 
kinetic energy of particles very effectively, addition of moderate interactions push the 
fractionally filled CI/QAH systems into the strongly correlated fractional topological 
phases. Some numerical works have systematically explored the interaction effects 
within the TFB models, and have firmly established the existence of FCI/FQAH 
states jsiEiniiniiis]. Recently, there are various proposals of material and cold-atom 
realization schemes for such FCI/FQAH systems [Ml EHl UHl HB EHl UHl 1201 Ell EB ESI ES] • 

For the continuum FQH system of two-dimensional (2D) electron gas in a strong 
magnetic held, Laughlin [25] used the lowest-Landau-level (LLL) single-particle wave 
functions with conserved angular momenta to construct the many-body wave function 
analytically for the FQH states on disk geometry. Then some other works pointed 
out the wave functions can also be constructed for continuous systems on different 
geometries like a sphere, a torus or a cylinder [26l EZ]- Later works found that 
the Laughlin wave function can be decomposed with the help of the Jack symmetric 
polynomial (Jacks) [281 EH], which naturally rehects the generalized Pauli principle 
(GPP) [HOI Ell EEl [33] in the Fock space for both the Abelian and the non-Abelian 
FQH states. Just as the Pauli exclusion principle corresponds to the Fermi statistics, 
the GPP corresponds to the fractional statistics of anyons. 

Recently, hnding the suitable trial wave functions of FGI/FQAH states becomes an 
outstanding problem. It is conjectured that there is a one-to-one mapping between IQH 
and GI/QAH states [IDl EB EHl EEl EZ]- other words, single-particle orbitals of an LL 
can be mapped to the ones in a TFB. Based on the one-dimensional (ID) maximally 
localized Wannier functions, Qi et al have used a mapping between the FQH and the 
FGI/FQAH states to construct generic wave functions of FGI/FQAH states on cylinder 
geometry with these ID Wannier functions [TOl EB EHl ESI EZ] • In terms of this mapping 
relationship, the Haldane pseudo-potential izg for these FGI/FQAH states can be 
constructed [38l EH] through a proper gauge choice for Wannier functions. An improved 
prescription has been adopted to construct variational wave functions of FGI/FQAH 
states on torus geometry by utilizing the gauge-hxed (non-maximally) localized Wannier 
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states HOI amis. From another aspect, conventional FQH states can also be obtained 
for 2D lattices analytically by using the conformal held theory ^31IH] 05] . In contrast 
to the above analytical or semi-analytical approaches, we here pursue a very direct yet 
effective purely numerical prescription to construct FCI/FQAH states on disk geometry, 
without any variational parameter or adjustable gauge freedom, but just utilizing the 
powerful GPP and the Jacks structure of FQH states, and the information from exact 
numerical single-particle orbitals of TFB models. 

In this paper, we exploit the single-particle wave functions of CI/QAH states (in 
the Kagome model and the Haldane [H [9] model) with TFB parameters on 

disk geometry, and explore the polynomial structure of the continuum Laughlin wave 
functions [25] to establish the many-body wave functions of the FCI/FQAH states. 
We hrstly construct the hnite-size lattices on disk geometry and diagonalize the non¬ 
interacting Hamiltonian matrix to obtain the single-particle CI/QAH wave functions, 
then use the GPP [26] and the Jacks [28l[28| to construct the FCI/FQAH wave functions 
for interacting spinless fermions. In order to test the feasibility of the many-body wave 
functions of FCI/FQAH states, we calculate the overlaps between these wave functions 
(WFs) and those from ED studies. Very high WF overlaps demonstrate the feasibility 
of our method even when the Hilbert space dimension is more than 3 x 10®. Besides, our 
method makes it possible to treat bigger lattices hlled with more fermions (the single¬ 
particle orbitals for lattice disks with 2000 sites can be easily constructed, and the Jacks 
can provide the basis states for systems with more than 15 particles). We also construct 
the edge-excitation wave-functions for FCI/FQAH states on disk by using the GPP 
and the information of single-particle orbitals. We demonstrate that the degeneracy 
sequences of the edge excitations of fermions on disk geometry are just the same as 
the bosonic systems [IT] 06] , in quite good agreement with the chiral Luttinger liquid 
theory 08] . 



Figure 1. (color online), (a) The kagome lattice and (b) the honeycomb lattice on 
disk. The kagome and honeycomb lattices both satisfy the Cq rotational symmetry on 
disk. Various disk sizes are indicated by the circles and the labelled numbers. 
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2. Single-particle CI/QAH orbitals on disk 


CI/QAH states are a class of topological insulators with non-trivial phases labeled by the 
Chern numbers |1]. Haldane has proposed the hrst CI/QAH model [3] on a honeycomb 
lattice in 1988. Recently, a series of CI/QAH models with TFBs were proposed [HliniiTj. 
We put the CI/QAH models, e.g. the Haldane and Kagome models, on the disks with 
the Cq rotational symmetry (Fig. [1]). Using tight-binding approximation and adjusting 
the hopping and phase parameters of CI/QAH model, we obtain the TFBs, i.e. the 
CI/QAH model with the specihc parameters which determines a flat energy band with 
a high flatness ratio as well as a non-zero Chern number. An additional trap potential 
is required on the finite-size disk to confine the FCI/FQAH droplet and generate a soft 
boundary, in order that the edge modes are able to propagate around these geometries 
then we can observe clear edge states. We now consider the Haldane and Kagome 
lattice models on disk as Ref. |16] . Here we choose the conventional harmonic trap with 
the form V = Vtrap Sr in which Urap is the potential strength (with the nearest- 

neighbor hopping t as the energy unit), |r| as the radius from the disk center (with the 
half lattice constant a/2 as the length unit) and rir as the fermion number operator 

[ISlIlTj. 


LU 



(b) 


(c) 




L 


L 


L 


Figure 2. (color online). Single-particle CI/QAH orbitals on disk with TFB 
parameters. Energy spectra versus the angular-momentum quantum number for (a) 
the honeycomb-lattice Haldane model with Ng = 96 sites, (b) the kagome-lattice model 
with Ns = 72 sites. The harmonic trap potential parameter is Urap = 0.005 for both 
models. The ground-state angular momentum quantum number is 5. (c) Low-energy 
orbitals of the kagome-lattice model. The numbers labelled under the energy levels 
stand for the angular-momentum quantum numbers (without modulo operation) and 
the orders of lowest-energy states, and the red balls represent the fermions occupying 
in these single-particle orbitals. 


We hrst notice that the Hamiltonian of the Haldane honeycomb (HC) lattice model, 
loaded with interacting spinless fermions 0116], can be written as 


^HC = -t' 

((rr'» 


/J,/rexp (t^r'r) + H.C. 
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and the kagome (KG) lattice model, also loaded with interacting spinless fermions, is 
written as |16] . 


^KG = - ^ E [/r'/r eXp (#r'r) + H. 

(rr') 

E [/r'/r + H.C 

({rr'» 


C. 


+ ViE nrUr' + V2 ^ nrUr' 

(rr'> {(rr'}) 


( 2 ) 


where /r(/r) annihilates (creates) a spinless fermion at site r; (...), ((...)) and (((.. .))) 
denote the nearest-neighbor (NN), the next-nearest-neighbor (NNN), and the next-next- 
nearest-neighbor (NNNN) pairs of sites. The phases 0r'r = ±0 which break the time 
reversal symmetry [3], and Vi, V 2 are the NN and NNN repnlsive potentials. In order to 
get a TFB, we choose t = 1, t' = 0.60, t" = —0.58, 0 = 0.47r [9] for the Haldane model 
and t = 1, t' = —0.19, 0 = 0.227r |46] for the kagome model. 

It is easy to obtain the single-particle eigen-energies and eigen-states of the Haldane 
and Kagome models on hnite-size disks by constrncting the lattice with the given 
hopping parameters and diagonalizing the non-interacting Hamiltonian matrix. In this 
paper, we hrst choose the honeycomb lattice with 96 sites and the kagome lattice with 
72 sites on disk and add a trap potential to restrain the movement of particles. The 
angnlar momentnm is obviously conserved in these hnite-size disk systems and these 
single-particle orbitals are shown in the sectors labelled by angular momentum quantum 
numbers [Fig. [2] (a) and (b)]. The angular momentum quantum number L varies 
from 0 to 5 because of the Cq rotational symmetry. With a harmonic trap potential 
Ktrap = 0.005, the angular momentum of GS energy is 5 and the others increase with 
L = 4 -|- i for the low energy states where i denotes the ith lowest energy states [Fig. [2] 
(c)]. By the way, we show that the GPP commands the particles occupying in these 
orbitals. For example, considering a 1/3 FGI/FQAH state, there are 7 particles hlled 
in the specihc orbitals. Based on the GPP for the 1/3 state, there is no more than 
one particle in three consecutive orbitals. In order to guarantee this system with the 
lowest energy (GS), 7 particles £11 in the 19 orbitals as shown in Fig. |2](c). At the same 
time, the density profile of the single-particle orbitals on lattice sites with a harmonic 
trap potential Vtrap = 0.005 can be obtained along the radius of the disk as shown in 
Fig. |3l The single-particle probabilities in the real space also satisfy the Cq rotational 
symmetry (shown in Fig. |3]) because of the lattice structure of the kagomGlattice disk. 
With these single-particle orbitals, the next step is to generate the many-body wave 
functions for the FGI/FQAH states. 
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Figure 3. (color online). The single-particle orbitals on the kagome-lattice disk 
with TFB parameters. The probability density of the single-particle orbitals on a 
kagome disk with Ns = 72 sites. The harmonic trap potential is Vtrap = 0.005, (a- 
b) ground-state density profile, (c-d) the first-excited-state density profile, (e-f) the 
second-excited-state density profile. 


3. Jack polynomial and FCI/FQAH wave fnnctions 


The Laughlin wave function is an analytic function to describe the continuum FQH 
states of 2D electron gas in a magnetic field. The single-particle orbital on disk 
is = (27rm!2™')“^/^z”* exp(—|zp/4) with angular-momentum quantum number 

m = 0,1, 2, 3,... (i.e. the angular momentum is = mh). The z/ = 1/3 Laughlin wave 
function on an inhnite-size disk [25] is 


i<j 


Ek^r/4/i 


(3) 


where z=x+iy and Ib is the magnetic length. We drop the non-universal Gaussian factor 
exp(— ^ and investigate the universal polynomial structures of = 

n(^* - 

i<j 

As a tutorial example, we consider the 1/3-Laughlin state with just three fermions, 
and expand the wave-function as follows: 

^2, ^ 3 ) = (D - Z2fizi - Z3f{z2 - Zsf 

= { + 1){ZI^Z2^Z3^ - Zi^Z2^Z3^ - Zi^Z2^Zs^ + Zi^Z2^Z3^ + Zi^Z2^Z3^ - Zl°Z2^^3®) 

+ (~3)(Zi®Z2^^3 ~ Z^Z2Z‘^ ZxZ’i’Z'^ Z\ Z2Z‘^ — Z\Z2^Z'^ — Zi^Z2^Z3) 

+ (-3)(zi^Z 2^Z3° - Zi^Z2°Z3^ - Zi^Z2^Z3° + Zi^Z2°Z3^ + Zi°Z2^Z3^ - Zi°Z2^Z3^) 

-h (-h6)(Zi^Z2^Z3 - Zi^Z2Z3^ - Zi^Z2^Z3 + “ ZiZ2^Z3^ + ^i^^2^3^) 

+ (-15)(z3Wz 3^ - Z3%W - ziWz3^ + + D^;^2S3" - D^^2=*;^3"). 

This expansion form of the Laughlin wave function can be classified by the exponents 
of anti-symmetric polynomials. For example, we can write the first term (zi^Z 2 ^Z 3 ^ — 
zi^Z 2 °Z 3 ^ - Z 1 ^Z 2 ^Z 3 ° + zi^Z 2 °Z 3 ^ + zi°Z 2 ^Z 3 ^ - zi°Z 2 ^Z 3 ^) as d'[ 6 , 3 , 0 ]- And [ 6 , 3 , 0 ] denotes 
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a configuration for three particles filled in three orbitals with the angular-momentum 
quantum number m = 0, 3 and 6 respectively. One can easily hnd that the 'h[ 6 , 3 , 0 ] 
is just an (unnormalized) Slater determinant. So the fit[ 6 , 3 , 0 ] can be construed as a 
fermionic many-particle state in which the three particles fulhll the Pauli exclusion 
principle. The three-fermion Laughlin state (shown in the above) can be written as 
^ 2 , ^ 3 ) = (+1)'I^[6,3,o] + (“3)T[6,2,i] + (—3)T[5^4^o] + (+6)T[5^3^i]-)-(— 15)T[4^3_2]- That 
is to say, we can decompose this three-fermion Laughlin state into hve (unnormalized) 
Slater determinants. However, it is very difficult to expand the Laughlin states with 
more than hve particles by such a brute-force method. 

The Laughlin state can also be decomposed into symmetric monomials (for bosons) 
or anti-symmetric Slater determinants (for fermions) with the aid of the Jacks [2811^ . 
The decomposition coefficients, e.g. -]-l, —3, —3,-f6, —15 in the above three-fermion 
Laughlin state, can also be generated very effectively by the Jacks. Such an approach 
even works for more than 15 particles. To be explicit, we decompose the fermionic 
FQH states (e.g. with the hlling factor v = 1/3) into Slater determinants by the Jack- 
polynomial method [281129]: 


= E (4) 

fi<X 

where A is the root partition which can index the basis state like X=[k0'^~^k0'^~^...k] 
which shows the sequence of occupation numbers of fermions in single-particle orbitals 
labelled by angular-momentum quantum numbers, denotes a Slater determinant 
and n stores the angular-momentum quantum numbers of each partition. For a particle 
hlling number z/ = there is a Jacks parameter a = —The expansion coefficient 


bxfj, in the Eq. 


b\^, = 


can be obtained by this recurrence relation 

nP'irA- n^(\] ^ (hi “ hi)&Ae(-l) 

P\ [^) Pfi 0 ; /i<e<x 


^sw 


(5) 


where Px{a) = E* K[Xi 2i{l - 1/a)]. Partitions 6=[pi, p 2 , + s, ~ s, ...,Pn] 

and p=[fii, p 2 , ■■■, Pi, ■■■, Pj, ■■■, Pn] are squeezed from the root partition A. Nsw is the 
number of swaps required to order the squeezed partition. We set bxx = 1 for simplicity. 

In the following, we will show the procedure to obtain the decomposition coefficients 
by squeezing the partitions successively. Consider the previous 1/3 Laughlin state with 
3 fermions as an example, the Jacks parameter is a = —2, and the root partition 
is A =[1001001] with the angular-momentum partition [6,3,0] which satishes the GPP. 
And three partitions Ai=[1000110] (with angular-momentum quantum numbers [5,4,0]), 
A 2 =[ 0110001 ] (with angular-momentum quantum numbers [ 6 , 2 , 1 ]), A 3 =[ 0101010 ] (with 
angular-momentum quantum numbers [5,3,1] )are squeezed from the root partition, but 
preserve the original total angular-momentum quantum number. The partition A 3 can 
be also squeezed from Ai or A 2 . The partition A 4 =[ 0011100 ] (with angular-momentum 
quantum numbers [4,3,2]) is squeezed from A 3 . The decomposition coefficients bx^ can 
be obtained recursively by the recurrence relation Eq. (I5|), are just the same coefficients 
J-l, —3, —3, -)- 6 , —15 which have been obtained by the previous brute-force method. And 



the 1/3 Laughlin state can be written by the Jacks in the second-quantization form 
as = JillOOlOOl) - J2IIOOOIIO) - J3IOIIOOOI) + J4IOIOIOIO) - J5IOOIIIOO) 

where J* is related to the decomposition coefficient -|-1, —3, —3, -|-6, —15 while now being 
normalized. The |1001...) corresponds to the normalized Slater-determinant state for 
fermions. Since there are hve different partitions (i.e. hve Slater-determinant states), 
and thus the squeezed Hilbert-space dimension Ujacks = 5 for this 1/3 Laughlin state 
with three fermions. 

Now consider the mapping between an LL and a TFB orbital on disk geometry, 
and we also assume that the many-body states of FCI/FQAH have similar polynomial 
structure as the continuum FQH state on disk. We can construct the many-body 
wave functions with the Jacks structure of the Laughlin FQH states while using the 
TFB single-particle orbitals (instead of the original LL orbitals). We have mentioned 
before that, the single-particle state in an LL is of the form 0m(^i) ~ The above 
Slater-determinant states, e.g. 11001001), with an angular-momentum partition [6,3,0], 
can be written as an unnormalized form, 11001001) rsj j-Q 2 Qj [ 06 (^ 1 ) 03 ( 2 ^ 2 ) 00 ( 2 : 3 ) - 

00(^1)00(2:2)03(2:3) - 03(2:1)06(2:2)00(2:3) + 03(2:1)00(^2)06(2:3) + 00(2:1)06(^2)03(2:3) - 

00 (^ 1 ) 03 ( 2 : 2 ) 00 ( 2 : 3 )]. Now based upon the one-to-one correspondence between LL and 
TFB orbitals (with the increasing sequence of angular-momentum quantum numbers), 
we can replace (j)m{zi) in the Slater determinants by the numerically obtained CI/QAH 
single-particle orbitals with the TFB parameters [Fig. [21(c)], then the Jack polynomials 
of the conventional FQH states can be translated into the many-body FCI/FQAH lattice 
states. Such a numerically obtained FCI/FQAH lattice state is labelled as Tjacks since 
it carries the same Jacks structure as the conventional FQH states. Just to keep 
things straight, there are 7 particles hlled in the TFB single-particle orbitals in Fig. 
He ) and we begin with this conhguration for the FCI/FQAH state with the minimum 
angular momentum for 1/3 Laughlin hlled based on the GPP. This conhguration can be 
recorded as 11001001001001001001) and it is the root conhguration that the others can 
be generated from it with the hxed same angular momentum using the above mentioned 
squeezing rules. The edge excitations of the Laughlin state on disk can be generated by 
the symmetric polynomials |18]. The root conhgurations for the hrst excited 

states can be written as 110010010010010010001). Based on the above procedures, the 
many-body lattice wave functions for the z/ = 1/3 FCI/FQAH states can be constructed 
numerically. 

The single-particle density operator is dehned as: 

N 

p(r) = ^J(r-r0 (6) 

i=l 

where r* denotes the position in a 2D coordinate, N is the particle number, and 
6{r — Vi) is a 2D delta-function. For a given wave function |\k(ri, r 2 ,..., r^r)), the 
density function is dehned as 

/ ^ ^ (T(ri, r2,..., rjv)|p(r)|T(ri, r2,..., rjy)) 

^ ^ (T(ri,r2,... ,rAr)|T(ri,r2,.. 


( 7 ) 
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Figure 4. (color online), (a) The real-space density profile of the continuum 1/3 
Laughlin state on disk. With the increasing of particles, p(0) —>■ 1/3. (b) The real- 
space density profile of the constructed 1/3 FCI/FQAH states dtjacks on the 180-site 
kagome-lattice disk. 


It is a difficult task to calculate r 2 ,..., rAr)|/5(r)|'k(ri, r 2 ,..., rAr)) for a 

many-particle system analytically. The Jacks provides a shortcut to obtain the 
occupation density. With the expansion coefficient b\^ of every partition, it’s easy to 
obtain the occupation probability pm in a specific orbital with the angular momentum 
quantum number m on disk or sphere geometries. The real-space density p(r) of the 
Laughlin state on disk can be written with the aid of the Jacks: 

d(^) = E exp(-rV2) Pm (8) 

where pm denotes the occupation density in the mth orbital and the factor 
exp(_^ 2 / 2 ) jg related to the single-particle orbitals of the 2D electron gas in a 
strong magnetic field. The real-space density profile of the 1/3 Laughlin GS on disk is 
shown in Fig. |l](a), with the increase of the particle number, p(0) —)■ Inspired by this, 

the real-space density profile of the FCI/FQAH states on disk can be obtained [shown 
in Fig. Ill[b)]. Here we choose the 180-site kagomClattice TFB model with the | hlling 
and a harmonic trap potential Vtrap = 0.005 along the radial direction. The density 
profile of the FCI/FQAH states exhibit similar structures as the Laughlin states,and 
also similar variation tendency when the number of fermions is increased. We should 
note that the particle density in the FCI bulk (i.e. the average number of fermions 
per lattice site) is around 0.11 as shown in Fig. |4l[b), since the filling number 1/3 for 
the FCI systems is for the filling of fermions in orbitals but not in lattice sites. The 
kagome lattice has three energy bands (due to a three-site unit cell), fermions here only 
effectively occupy the single-particle orbitals of the lowest energy band (i.e. the TFB), 
and thus the particle density in the bulk of the 1/3 FCI states should be close to the 
value (1/3) x (1/3) = 1/9. 
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Figure 5. (color online). Comparison of wave functions constructed by the Jacks with 
that of the ED studies, e.g. O = |(4'ED|'I'jacks)|- In these figures, the red solid circle 
represents the overlap between the wave functions of 1/3 FCI/FQAH states constructed 
by the Jacks with the ED results for different repulsion interactions Vi. (a)-(c) The 
GS of the 42-site kagome-lattice disk filled with 2,3 and 4 fermions, (d) The GS of the 
54-site kagome-lattice disk filled with 5 fermions, (e)-(f) The first excited states of the 
54-site kagome-lattice disk filled with 3 and 4 fermions. The harmonic trap potential 
is set as Ttrap = 0.005 as before. 


4. Comparison of the wave functions 

It is almost impossible to write down an exact wave function of FCI/FQAH states 
analytically like the Laughlin wave functions because of the discontinuity of the CI/QAH 
lattice systems in contrast to the continuity of 2D electron gas. Though we cannot 
describe FCI/FQAH states by an analytical expression exactly, the wave functions of 
the FCI/FQAH states can be written as ^ by the Jacks [28] 

i 

where Ji are normalized Jacks coefficients of the corresponding FQH state and 
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Table 1. The overlaps, O = |('I'ED|'l'jacks)|, between the 1/3 FCI/FQAH states 
constructed from the Jacks and that of the ED calculations for the kagome-disk. 
The overlaps for the ground states (Ogs) and the first excited states (Oes) with 
various fermion numbers (-/V/’s) and lattices sites (A'l’s) are listed. The Hilbert- 
space dimensions of the ED calculations (Ved) and that of the Jacks ("ZJjacks.GS for 
the ground-state Jacks, Ujacks.ES for the first-excited-state Jacks) are also listed for 
comparison. 


Nf 

Nl 

2^ed,gs 

^Jacks,GS 

Dgs 

Nl 

TJed.es 

7Jjacks.ES 

Oes 

2 

42 

861 

2 

0.998 

42 

861 

2 

0.999 

3 

42 

11480 

5 

0.995 

42 

11480 

6 

0.996 

4 

42 

111930 

16 

0.984 

42 

111930 

20 

0.980 

5 

54 

3162510 

59 

0.971 

66 

8936928 

77 

0.903 


denotes the Slater determinants. Instead of the IQH state which obeys the Pauli 
exclusion principle, the particles in the FQH system obey the GPP. From the ED 
studies on the torus geometry, we know that it is necessary to add the interaction 
potential Vi on the NN bonds to obtain the 1/3 FQH phase. In order to estimate the 
rationality of the many-body wave functions constructed by using the Jacks, we can 
compare it with the ED results of wave functions in the real space, e.g. calculating 
O = |{J/ED|d'jacks)| for FCI/FQAH states. We hrst consider the FCI GS wave function 
of the kagome lattice and the results of the WF overlaps are shown in Fig. |5](a)-(d). For 
the large enough interaction potential, we can observe a high overlap. Though the Jacks 
provides a much fewer squeezed Hilbert-space dimension than the ED dimension (show 
in Table. [I]), the values of the overlaps for the kagome lattice are quite high. It is worth 
emphasizing that the maximum value of overlaps between the ED WFs and the Jacks 
WFs is 0.998, 0.996, 0.984, 0.971 for the systems with Nf = 2, 3, 4 and 5 fermions. We 
hnd that a good WF overlap might be also observed even without the trap potential. 
When the trap potential parameter Idrap varies from 0.0 to 0.0065 for the above cases, 
the overlap becomes increasingly higher. However when the trap exceeds 0.0065, the 
overlap decreases, and the best WF overlap needs an intermediate trap potential. 

With the help of the GPP, the excited states can also be constructed. It is very 
easy to obtain the hrst excited states in the FCI/FQAH phases like that of the GS. Here 
we show the overlap (Fig. [5] (e) and (f)) between the exact excited state wave function 
of kagomClattice FCI/FQAH state and the constructed Jacks states. The maximum 
value of the WF overlaps is 0.999, 0.996, 0.980, 0.903 for 2-, 3-, 4- and 5-fermion systems 
respectively (show in Table. H]). By the way, the overlaps of the kagomClattice excited 
states are also very robust when varying the repulsion interaction Vi. 
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5. Edge excitations on disk 

Edge excitations can be exploited to extract the crucial information of topological 
order of some continuum FQH liquids [HI 091 00]. In a previous work |16|, we 
investigated the edge excitations of the bosonic FCI/FQAH phases [9] on a disk geometry- 
using systematic numerical exact diagonalization (ED) studies and obtained the edge 
excitation spectra which is quite coincident with the chiral Luttinger liquid theory [H] . 
ED provides an accurate tool to deal with FCI/FQAH edge excitations, however it is 
difficult to study a large lattice system with more than a hundred sites and/or with 
more than several particles due to the Hilbert space limitation. The edge excitations of 
the Laughlin state on disk can be generated by the symmetric polynomials 
08] . and recently the Jacks has been employed to construct the edge states for the 
continuous FQH systems EH. 
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Figure 6. (color online). Edge excitations of the FCI/FQAH states on the 90-site 
kagome-lattice disk filled with various fermions, Hrap = 0.005. The labeled numbers 
show the quasi-degeneracy of the low-energy excitations, (a)-(d) show us the quasi¬ 
degenerate sequences for 2, 3, 4 and 5 fermions respectively 


In the above section, we have compared our FCI/FQAH wave functions with the 
ED results and the high overlaps verihed the feasibility of constructing many-body wave 
functions by the GPP and the Jacks. We will here try a direct approach to estimate the 
relative many-body energies for the FCI/FQAH states by the Jacks and the GPP. For 
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the 1/3 Laughlin states, particles in the single-particle orbitals fulhll the GPP. That is 
to say, there is no more than one particle to occupy any three adjacent orbitals. As we 
all know, for the identical fermion systems, the wave functions can be generated by the 
Slater determinant and the approximate energy of systems can be viewed as summing all 
energies of every single-particle orbital. For the FQH states, the many-body states are 
denoted by the Jacks instead of the Slater determinant. But we can expand the Jacks 
based on the Slater determinants. The approximate energies of the Laughlin systems 
can be calculated with the aid of the GPP and the Jacks: 

A'l = X! (^) 

i 

where pi is the probability of every basis related to the coefficient of Jacks and Ei is the 
energy of every conhguration. Now, we consider the kagomGlattice TFB model with the 
typical harmonic trap potential Idrap = 0.005 along the radial direction on the disk |16]. 
We put some particles in the single-particle orbitals (shown in Fig. [21(c)). As a quick 
case study, we put 3 fermions in the TFB and the behavior of particles obeys the GPP. 
The GS energy is composed by .F|iooiooi>,-E|ioooiio),-E|oiioooi>,-E|oioioio),-E|ooiiioo), where 
F'liooiooi) denotes the energy of the state 11001001). i?|iooiooi) is the sum of single-particle 
eigen-energies for fermions occupying the hrst, forth and seventh orbitals in Fig. He). 
Just as the excited-state WFs, the approximate excitation energies can be generated by 
the Jacks and the symmetric polynomials |18]. By this, we can obtain the approximate 
many-body energies of any FGI/FQAH state. Some representative results are shown 
in Fig. [6l Each state can be classihed with the angular momentum quantum number 
L = 0,1, 2, 3, ...mod 6 since the kagome lattice has the Cq rotational symmetry. The 
quasi-degeneracies of the edge excitations in fermionic FGI/FQAH systems are the same 
as that in the bosonic FGI/FQAH systems [47l l46] . 

Estimation of the energies using the Jacks above can be viewed as estimating the 
kinetic energy term of the Hamiltonian, and the repulsive interaction Vi term has not 
been taken into account explicitly. For a comparison with the ED results with hnite 
Vi values, we resort to an extrapolation scheme to effectively remove the repulsive 
interaction term by extrapolating the ED data in the FGI/FQAH region (with good 
WF overlaps and beyond some critical Vi values as shown in Fig. |5]) to the virtual 
Vi = 0 limit. An exponential htting can be chosen to ht the variations of the ED 
energies with different repulsive potentials. For 42 sites with two fermions, the absolute 
error of the ground-state energy is 0.000003978. For 42 sites with three fermions, the 
absolute error is 0.005233 for the ground-state energy, and is 0.004966 for the hrst- 
excited-state energy. For 54 sites with three fermions, the absolute error is 0.006004 
for the ground-state energy, and is 0.006129 for the hrst-excited-state energy. For 54 
sites with hve fermions the absolute error of ground-state energy is 0.01024. All of the 
relative errors of ground-state energies are less than 0.1% for these systems. 

The GPP and the Jacks provide us an alternate method to obtain the edge 
excitations of the FGI/FQAH states on disk geometry without resorting to the ED 
method. We can diagonalize a very large lattice to get the all single-particle eigen- 
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Table 2. The quasi-degeneracy sequences of FCI/FQAH states on disk with various 
fermion numbers. This degeneracy is obtained by the GPP and the Jacks. Here, dn 
marks the sectors of the degenerate energy levels, Nf is the number of fermions. 
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energies. There is a recurrence relation of edge-excitation quasi-degeneracy sequence for 
an infinite lattice with the finite-particle counting shown in Table. HI By utilizing this 
table, the degeneracy sequence can be read easily. We also can estimate the number 
of the particles and the size of the lattice which can generate a specific degeneracy 
sequence. By the way, the degeneracy is not related to the types of particles, and it 
is only controlled by the GPP. Maybe for a hnite-size lattice system hlled with more 
particles, different degeneracy sectors may cross their boundaries then the degeneracy 
sequence will be broken. 

There is some obvious advantage for utilizing the GPP and the Jacks to obtain 
the quasi-degeneracy sequence. Gomparing to the costly ED method, this computation 
method not only is simple but also can reach very large lattice size. In Fig. [71 we 
show the degeneracy sequence for some over 200-site lattice systems with 8 fermions 
hlled on disk which the ED method is not able to compute. By the way, the Haldane 
honeycomb-lattice model and the kagome-lattice model own the Cq rotational symmetry 
and the checkerboard-lattice model owns the C 4 rotational symmetry on disk geometry, 
the period of the angular momentum quantum number is 6,6,4 for the Haldane model, 
the kagome model and the checkerboard model, respectively. Though they own different 
rotational symmetries, their degeneracy sequences are identical. 
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Figure 7. (color online). Edge excitation in three large-size FCI/FQAH disks 
with 8 fermions. The labeled numbers show the quasi-degeneracy of the low-energy 
excitations, (a) The honeycomb-lattice Haldane model with 228 sites . (b) The 

kagome-lattice model with 228 sites . (c) The checkerboard-lattice model with 216 
sites. The same trap potential Ttrap = 0.005 has been chosen for all three models. 


6. Summary and discussion 

We use the numerical single-particle CI/QAH wave functions to generate the many- 
body wave functions for fermionic FCI/FQAH states on disk geometry. The single¬ 
particle CI/QAH eigen-states eigen-energies can be obtained by constructing a finite 
TFB lattice system on disk with harmonic trap potentials and diagonalizing the single¬ 
particle Hamiltonian numerically. We then utilize the GPP and the Jacks to generate 
the numerical many-body wave functions. Such approach yields a very direct yet 
effective purely numerical prescription to construct FCI/FQAH states on disk geometry, 
without any variational parameter or adjustable gauge freedom. In order to check the 
reasonability of our method, we make a WF overlap between the Jacks and the ED 
results. High values of the WF overlaps and the robustness against interaction conhrm 
the rationality. The edge excitations can also be approximately obtained with the single¬ 
particle orbitals and the GPP. The quasi-degeneracy sequences of fermionic FCI/FQAH 
systems reproduce the prediction of the chiral Luttinger liquid theory, and they are the 
same as previously studied the bosonic FCI/FQAH systems, complementing the exact 
diagonalization results with larger lattice sizes and more particles. 

We would like to mention some possible future research directions following along 
the present work. It is possible to further explore or utilize the FCI/FQAH wave 
functions at hand now, e.g., compare them with the lattice Laughlin states from 
the conformal held theory |l3l [HJ HS] and the Gutzwiller-projected parton wave 
functions [521 ESI El [55], and extracting the quasi-particle fractional statistics via the 
modular matrices [56l EH El EH]- It would be also very interesting to extend the present 
approach to the non-Abelian FGI/FQAH states [601 El [62], the FGI/FQAH states in 
high-Ghern-number TFBs [631E1E5], and also the hierarchy FGI/FQAH states [661 [67] . 
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